S'-bases as a tool to solve reduction problems for Feynman integrals 
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We suggest a mathematical definition of the notion of master integrals and present a brief review of algorithmic 
methods to solve reduction problems for Feynman integrals based on integration by parts relations. In particular, 
we discuss a recently suggested reduction algorithm which uses Grobner bases. New results obtained with its help 
for a family of three-loop Feynman integrals are outlined. 



1. Introduction 

In the framework of perturbation theory quan- 
tum-theoretical amplitudes are written in terms 
of Feynman integrals that are constructed accord- 
ing to Feynman rules. After a tensor reduction 
based on projectors is performed, each Feynman 
graph generates various scalar Feynman integrals 
with the same structure of the integrand and with 
various powers of propagators (indices): 

F{ai,...,an) = y'-'y-^ai ^H;r- (1) 

Here hi, i — 1, . . . , /i, are loop momenta and the 
denominators Er are either quadratic or linear 
with respect to the loop momenta pi — ki, i — 
1 . . . . , /i or to the independent external momenta 
Ph+i = qi,---,Ph+N = Qn of the graph. Ir- 
reducible polynomials in the numerator can be 
represented as denominators raised to negative 
powers. Usual prescriptions k'^ = k'^ + iO, etc. 
are implied. The dimensional regularization with 
d = 4 — 2£ is assumed. 

In today's perturbative calculations, when one 
needs to evaluate thousands or millions of Feyn- 
man integrals , a well-known traditional strat- 
egy is to apply integration by parts (IBP) rela- 
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tions T" and some other relations in order to find 
an algorithm that expresses any given Feynman 
integral as a linear combination of some master 
integrals. Then the whole problem of evaluation 
of Feynman integrals of a given family is decom- 
posed into two parts: constructing a reduction 
procedure and evaluating master integrals. 

It turns out, however, that there is no common 
definition of master integrals. Practically, after 
solving the reduction problem for a given family, 
we know that these are master integrals because 
we see them. There are examples where different 
authors consider different master integrals for the 
same class of Feynman integrals. In the next sec- 
tion, we are going to suggest a definition of this 
notion. In Sections 3 and 4, we give a short re- 
view of algorithmic methods to solve recursion 
problems in the framework of this definition. In 
particular, we characterize in Section 4 our algo- 
rithmic method |2I3| which uses so-called s-bases 
and outline new results obtained with its help for 
a family of three-loop Feynman integrals. 

2. The notion of master integrals 

To define the notion of master integrals we need 
some definitions. Feynman integrals Q can be 
considered as elements of the field of functions T 
of n integer arguments oi, 02, . . . , a„. This is an 
infinitely dimensional linear space. The simplest 
basis of this space is the set of elements i?ai,...,a„, 

where i?ai,...,a„(ai, . . . ,a^) = Sa^y^ . . - Sa^^a'^ ■ 
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One usually considers some set of relations be- 
tween Feynman integrals. Formally, they are el- 
ements of the adjoint vector space T* ^ i.e. the 
linear functionals on T , so that, for any r G T* ^ 
there is a corresponding value (r, /) for any given 
/ G . The simplest basis of this space is the set 
of elements _ff*^ which are defined as follows: 

{Hau-,a„J) = /(ai,...,a„) . 

There are several types of commonly used rela- 
tions. The first and the main type consists of the 
integration by parts relations P 

J--Jii^''^^'i;{p^mA-o. (2) 

i' \ 1=1 / 

After differentiating, the scalar products ki ■ kj 
and ki ■ qj are expressed linearly in terms of the 
factors Ei of the denominator, and one obtains 
the IBP relations in the following form: 

^aiF{ai +bis, . ■ ■ ,an + b^^n) = , (3) 

Now one can substitute all possible 
(fli, 02, ... , a„) on the left-hand sides of (PJ and 
obtain an infinite set of elements of T* . 

In addition to the IBP relations one often con- 
siders the so-called Lorentz-invariance identities 
0]. It is not yet clear whether they follow from 
IBP relations but they turn out to be useful. 

Another type of relations comes from the sym- 
metries of the given family of integrals. Typically 
they have the following form: 

i^(ai, . . . ,a„) = (-1)^ i^(a^(i),...,a^(„)), 

where di are fixed and are equal to either one 
or zero, and tt is a permutation. Again, after 
substituting all possible we obtain an infinite 
set of elements of T*. 

Then one takes into account boundary condi- 
tions, i.e. the conditions of the following form: 

F{ai, . . . , a„) = when a^j < 0, . . . Qi^ < (4) 

for some subset of indices. 

One more type of relations corresponds to par- 
ity conditions. For example, Feynman integrals 
can be zero if the sum of some subset of indices 
is odd and each index is nonpositive. 



After having fixed a set of relations we can 
generate by them an infinitely dimensional vec- 
tor subspace TZ C J-*. Now one considers the set 
of solutions of all those relations, that is the in- 
tersection of the kernels of all functionals r G TZ. 
This is a vector subspace of that will be de- 
noted with S. A Feynman integral considered as 
a function of the integer variables ai, . . . , a„ is an 
element of the space S for it satisfies the IBP re- 
lations and other relations mentioned above. For- 
mally, 5 = {/ G ^ : (r, /} = V r G 7^}. 

The dimension of 5 might be infinite. However, 
for all known classes of Feynman integrals, where 
the reduction problem has been solved, it appears 
to be finite. It looks like it can be a non-trivial 
mathematical problem to prove that this holds 
for all families of Feynman integrals. 

When talking about expressing one Feynman 
integral by another, it is usually assumed that 
we consider the consequences of relations TZ. Let 
us say that an integral F{ai, . . . , a„) can be ex- 
pressed via some other integrals i^(aj, . . . , a^), 
. . ., F{ai, . . . , a^) if there exists r GTZ such that 

(r, F) F(ai, . . . ,a„) 

Let us turn to the notion of irreducibility of 
Feynman integrals. Suppose we have two inte- 
grals F{ai, . . . , a„) and F{a[, . . . , a'^) that can be 
expressed one by another, for example, due to a 
symmetry of the diagram. Of course, it is rea- 
sonable to choose only one of them as a master 
integral. However there seems to be nothing nat- 
ural in this choice, for they are equivalent. So, 
even having fixed a set of relations, we do not 
have enough information to define master inte- 
grals. The only thing wc know that their number 
is equal to the dimension of S. 

Therefore, to define master (or, irreducible) in- 
tegrals, we need to choose a certain priority be- 
tween the points (oi, . . . , a„), formally, to intro- 
duce a complete ordering on them (that will be 
denoted with the symbol -< and named as lower). 
There are different ways to do that but at least 
it looks natural to have simpler integrals corre- 
sponding to the minimal elements in this order- 
ing. We shall introduce an ordering in two steps. 
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First of all, let us realize that the Feynman inte- 
grals are simpler, from the analytic point of view, 
if they have more non-positive indices. In fact, 
in numerous examples of solving of IBP relations 
by hand, the natural goal was to reduce indices 
to zero or negative values.'^ The big experience 
reflected in many papers has led to the natural 
idea to decompose the whole region of the integer 
indices which we call sectors.^ 

Let us consider the set V of elements 
{di,...,d„} called directions, where all di are 
equal to 1 or —1. For any given direction 
h' — {di,...,dn}, we consider the region cr^ = 
{(ai,...,a„) : (oj — l/2)di > 0} and call it a 
sector. In other words, in a given sector, the 
indices corresponding to ±1 are positive (non- 
positive) . Obviously the union of all sectors con- 
tains all integer points in the n-dimensional vec- 
tor space and the intersection of any two sec- 
tors is an empty set. We will say that a di- 
rection {di, . . . dn} is lower than {d'l, . . . djj} if 
^^1 l£ d[, . . . , dn ^ d'^ and they are not equal. The 
same is said about the corresponding sectors. 

It looks natural to define that F(ai, . . . , a„) >- 
F{a[, . . . , a^) if the sector of {a[, . . . , a^) is lower 
than the sector of (oi , . . . , a„) . Moreover it is also 
natural to assume that the 'corner point' ((o?i + 
l)/2,...,(d„ -I- l)/2) (those numbers are either 
ones or zeros) of the sector cr{di,...,d„} is lower 
than all other points of this sector. 

To define an ordering completely we introduce 
it in some way inside the sectors. After this, 
we can define what a master integral is. It is 
such an integral F{ai, . . . , a„) that there is no el- 
ement r G TZ acting on F according to relation 

10) such that all the points {a[, . . . ,a'jj are lower 
than (ai, . . . , a„). 

As it has been noted in 7 , one can prove that 
a given integral F{ai, . . . , an) is a master integral 
by constructing a function c € T such that (*) 
i) it satisfies all relations TZ and therefore c G 5; 

11) c(ai, . . . ,a„) = 1; 



^Moreover, the reduction procedure often stopped once 
one arrived at integrals which are already sufficiently sim- 
ple, in particular, when they can be expressed analytically 
in terms of gamma functions for general d. 
*This decomposition is also standard in the so-called La- 
porta's algorithm |S] which is briefly characterized below. 



iii) c{a'i, . . . , a'n) = for all lower (a'^, . . . , a^). 
To do this one can use the basic parametric rep- 
resentation jnj of the Baikov's method [Jj which 
is discussed below. 



3. Algorithmic methods of solving 
reduction problems 

Now we can give a short overview of the known 
algorithmic methods from the standpoint of our 
definitions. 

The Laporta's algorithm j2j works with sys- 
tems of equations for individual Feynman inte- 
grals. Let Tt be the set the subspace of gen- 
erated by Haj^^,,,^a„ where \ Ui \< t and TZt be 
the intersection of TZ with the subset of JF* gen- 
erated by H*^ where \ ai\<t. The limit of 
the difference between the dimensions of Tt and 
TZt when t tends to infinity is the dimension of 
S, so that there is a certain t such that TZt has 
'enough' relations to express any given integral 
i^(ai, . . . , a„) with \ ai \< t as a, linear combi- 
nation of master integrals with the use of those 
relations. So, the method is based on finding such 
a t and solving a huge system of linear equations. 

There is a public version jBj of the implemen- 
tation of this algorithm as well as a number of 
private versions. 

The basic tool of the Baikov's method |7I8| (see 
also Chapter 6 of |H]) is the representation 

J J Xi . . . Xn 

where P is constructed for a given family of in- 
tegrals according to some rules. This represen- 
tation is used to prove that an integral with 
Oil, ... , OLin = is a master integral. Indeed, such 
a function automatically satisfies two of the three 
points of (*) and it is left to check whether it does 
not vanish at (oi, . . . , a„). 

Suppose that we know that the integrals F{A') 
with A"" = {a\,. . . ,a^), A'' = (aj, . . . ,aj;) are 
master integrals for we have constructed the cor- 
responding solutions of this type C(ai,...,a„) satis- 
fying the conditions (*). These functions form a 
subbasis of the solution space S. Let us suppose 
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that F can be represented as 
F — ^ ^ ki Ci . 



(7) 



(This is always true if we know the complete set 
of master integrals.) Then one substitutes all Aj 
and solves an upper-triangle system of equations 

F(^j) = ^A:,C,(Aj)fori<j. 



so the coefficients ki are expressed in terms of 
the values F{Ai). Then the knowledge of both ki 
and d provides the possibility to calculate any 
given Feynman integral of this class. In fact it 
can happen that there is more than one master 
integral in a given sector, therefore one should 
consider also points different from sector corners. 

This method has not yet been analyzed from 
the mathematical point of view. After construct- 
ing a set of master integrals it is not clear whether 
this is a complete set. Still one can always per- 
form various crucial checks and then understand 
that the recursion problem was solved indeed. 

4. Using Grobner bases and s-bases 

Tarasov suggested the idea ,11 to apply 
Grobner bases ^Hl to solve reduction problems 
for Feynman integrals^. In his approach, IBP re- 
lations are reduced to differential equations. An 
attempt to use Grobner bases associated with the 
shift operators was made in |16| . 

In the recently introduced method |2l3j which 
uses so-called s-bases (Grobner-type bases), the 
ordering in a given sector should be chosen more 
strictly. Let us fix a sector cr{ci....,c„) and take 
Pi = (ci -I- l)/2. Then let us say that a (multi)- 
degree of a point (ai,...,a„) e o'(ci,...,c„) is 
(ci(ai -pi), . . . ,c„(a„ -Pn))- The degree of the 
corner of the sector (pi, . . . ,p„) is obviously equal 
to (0, ...,0). The degrees form the semi-group 
(with respect to (6i, . . . , 5„) + (6i, . . . , 5'„) 

^As an application of this approach, the solution of the 
reduction problem for two-loop self-energy diagrams with 
five general masses was obtained in 1121 . with an agree- 
ment with an earlier solution M'M . The solution of the 
reduction problem for massless two-loop off-shell vertex di- 
agrams (which was first obtained in fI3l within Laporta's 
algorithm) was reproduced in 1151 . 



{hi+h'^, . . . ,hn + h'^)). An ordering on the points 
is naturally extended to the ordering of degrees. 
We will say that an ordering on N" (denoted with 
the symbol )^) is proper if 

i) for any a £ N" not equal to (0, . . . 0) one has 

(0, ...0) 

ii) for any a,b,c G N" one has a ^ 6 if and only 
if a + c ^ b + c. 

Now we are allowed to take any ordering such 
that the corresponding ordering on degrees is 
proper for each sector. This corresponds to the 
basic idea of the method to introduce the algebra 
acting on J- and to work with operators and their 
degrees. Let us outline it briefiy. This algebra A 
is generated by shift operators Yi and multiplica- 
tion operators Ai that act via 

(Y^ ■ F)(ai, . . . ,a„) F(ai, . . . , + 1, . . . , a„) , 
{A^ ■ F)(ai, . . . ,a„) = aiF{ai, . . . ,a„) . 

The IBP relations are interpreted as elements of 
this algebra and generate an ideal I. There is an 
algorithm |2l3j designed to build a set of special 
bases of this ideal, allowing to reduce any integral 
to master integrals. To elucidate, any element 
of this ideal corresponds to an infinite number of 
relations r € TZ obtained by substituting different 
Qi. Obviously, 



F{ai, . . . ,a„) 



Y[iY;^p(-.-P.) .f] (pi,...,p„). 



and the exponents at Y^^ form exactly the de- 
gree of (oi, . . . , a„). Hence expressing an integral 
with lower integrals corresponds to the idea of 
the method to build the elements of T with low- 
est possible degrees. 

If one has constructed the s-bases for all sec- 
tors, the reduction algorithm can be applied. It 
expresses any integral as a linear combination of 
a certain number of integrals, that could not be 
reduced by the method. 

This algorithm can work for problems with 
many indices. In [T7j . it was applied to three- 
loop Feynman integrals 



F{ai, ... ,09) 



d'^kd'^ld'^r 



{-2v ■ kY^{-2v ■ lY'^ 
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{2v ■ 

Here ag < 0, v is the HQET quark veloc- 
ity. The integrals are symmetric with respect 
to (1 ^ 2,3 ^ 4,7 <^ 8) and vanish if one 
of the following sets is with non-positive indices: 
05,7, 05,8, 06,7, a6,8, 07,8, 03,4,6- The algorithm of 
|2I3| has revealed the following master integrals: 

h = 1,0,1, 1,1, 1,0,0), 

h = F(l, 1,1, 1,0,0,1, 1,0), 

/3 = F(l, 1,0,0,0,1, 1,1,0), 

h - F(0, 1,1, 0,1, 1,0, 1,0), 

h = F(-l, 1,1, 0,1,1, 0,1, 0) , 

h = F(0, 0,0, 1,1, 1,1, 0,0), 

h = F(0, 1,0, 0,0, 1,1, 1,0), 

It = F(0, 1,0,0, 1,1,1,0,0) , 

It = F(0, 2, 0,0, 1,1,1,0,0) , 

h = F(0, 0,0, 0,0, 1,1, 1,0). 

These are some examples of reduction: 
3d - 10 f 3(d-4) 



F(1,...,1,0) 



d-5 
3(ci-4) 



h - 



8{2d - 9) 
(rf-3)(3d-8) 



16(2d-9) 8(3d- 13)(3d- 11)^ 
3(d- 2)(3d- ll)(3d- 8) 



F(l,...,l,-1) = 



64(2d- 9)(2d- 7)(3rf- 13) 

9(d-4)(rf-2)(3rf-8) 
'64(2d- 9)(2d- 7)(3rf- 13)' 

3(3<i-8) - 
~32(2d- 9)(2rf- 7) ^ 

3(d-3)(3d- 11) 



16(d-5)(d-4)(2d-9) 
(d- 2)(2d- 7)(2d- 5) ^ 



8(d- 3)(2d- 9)(3d- 13) 

3(2d - 7)2(2d - 5)(3d - ll)(3d - 7) 



256(d-4)2(d-3)(2(i-9) 

Preliminary checks have shown that the present 
algorithm can work for large ri, e.g., for n = 12. 
There are various interesting practical and math- 
ematical problems connected with method. We 



hope to solve them in the future and improve our 
algorithm for more sophisticated calculations. 
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